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We study theoretically the mixed state properties of a strong uniaxially-anisotropic type II super¬ 
conductor with the Pauli paramagnetic effect, focusing on their behaviors when the magnetic field 
orientation is tilted from the conduction layer ab plane. On the basis of Eilenberger theory, we quan¬ 
titatively estimate significant contributions of the Pauli paramagnetic effects on a variety of physical 
observables, including transverse and longitudinal components of the flux line lattice form factors, 
magnetization curves, Sommerfeld coefficient, field distributions and magnetic torques. We apply 
these studies to Sr2Ru04 and quantitatively explain several seemingly curious behaviors, including 
the Hc 2 suppression for the ab plane direction, the larger anisotropy ratio and intensity found by the 
spin-flip small angle neutron scattering, and the first order transition observed recently in magneto¬ 
caloric, specific heat and magnetization measurements in a coherent and consistent manner. Those 
lead us to conclude that Sr2Ru04 is either a spin-singlet or a spin-triplet pairing with the d-vector 
components in the ab plane. 

PACS numbers: 74.25.Uv, 74.70.Pq, 74.25.Ha, 61.05.fg 


I. INTRODUCTION 

Sr2Ru04 is well known to be a prime candidate of a 
chiral p-wave superconductor m- The crystal struc¬ 
ture is same as in La2Cu04: a mother compound of 
high Tc superconductors. The normal state properties 
of Sr2Ru04 are characterized by quite a standard Lan¬ 
dau Fermi liquid picture with a moderate mass renor¬ 
malization [T] in stark contrast with the high Tc cuprates 
which are strange metals in every respect. Yet both have 
a strong two dimensional metallic conduction associated 
with anisotropic layered structure. In this sense Sr2Ru04 
has a firm foundation, out of which the superconducting 
state develops at Tc=1.5K. Thus we can safely employ 
a reliable theoretical framework such as Eilenberger the¬ 
ory that assumes a normal Fermi liquid for describing the 
superconducting properties under an applied field. 

Recently, the research front of Sr2Ru04 has been 
greatly advanced: (1) The small angle neutron scattering 
(SANS) experiment [4] shows that the anisotropy ratio 
of the vortex lattice amounts to Evl ^ 60 for the field 
orientation B parallel to the ab plane. This is at odds 
with the Hc 2 anisotropy ratio T ^^2 = ^c2,ab/^c2,c = 20, 
where iLc2,ab (^c2,c) is the upper critical field Hc 2 for 
B II ab (B II c), because in usual single-band supercon¬ 
ductors Tvl ^ i^ expected [5]. 

(2) The magneto-caloric effect [6], the specific heat [7] 
and magnetization experiments [8] detect the first order 
transition at iLc2,ab iii low temperatures, which is similar 
to superconductors with strong Pauli paramagnetic effect 
(PPE), such as in CeCoIns [OUT^. 

We note that the three experiments [6]-[8] mentioned 
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above are mutually quite consistent with each others, 
since a certain amount of quasi-particles in the super¬ 
conducting state are responsible for exhibiting the first 
order transition’s jumps at iLc2,ab in those thermody¬ 
namic quantities. There the same quasi-particles man¬ 
ifest themselves in each observable. This means that 
viewing from the normal side above iEc2,ab the spin sus¬ 
ceptibility Xspin must decrease in the superconducting 
state. However, this expectation is in conflict with the 
existing Knight shift experiments by NMR [isHIl] and 
polarized neutron scattering measurements [16]. There 
are no triplet pairing theories proposed so far which are 
able to predict the first order Hc2,ab transition, including 
works by one of the present authors [T7H38] . 

There are some other experimental reports suggest¬ 
ing the spin-triplet chiral p-wave superconductivity in 
Sr2Ru04. For example, the observation of half-quantized 
fluxoids [39] , which requires multiple order parameter for 
the pairing function with both spin and orbital degrees 
of freedom active, implying the possibility of the spin- 
triplet pairing. The chiral domain formations and the 
time reversal symmetry breaking are suggested by vari¬ 
ous experimental methods [4QH42] . However, a scanning 
Hall probe experiment [43] fails to detect the edge current 
expected for the chiral superconductors. In the experi¬ 
ments, the estimated domain sizes for each sample used 
in those experiments are strangely widely different from 
Ipm to 1mm (see Ref. [44] for detailed critical exami¬ 
nations on this point). Therefore, in the present status 
of understanding the mechanism of superconductivity in 
Sr2Ru04, the above-mentioned experimental results are 
mutually contradicted. 

The purpose of this paper is to find a clue to resolve 
the contradictions, by describing the mixed state prop¬ 
erties for a uniaxial strong anisotropic type H supercon¬ 
ductor with PPE in the clean limit and a single band on 
the basis of quasi-classical Eilenberger theory. Then we 
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critically examine several experiments done recently on 
Sr 2 Ru 04 and interpret the implications of those experi¬ 
ments in the viewpoint of the PPE. It is shown that the 
results are maximally consistent with the experimental 
data, and stimulate future theoretical and experimental 
studies to further understand the mechanism of the ex¬ 
otic superconductivity. 

The Eilenberger theory is applicable for superconduc¬ 
tors with /cfC ^ 1- For Sr 2 Ru 04 this condition is well 
satisfied because the coherence length ^ ^ 30nm and the 
inverse of the Eermi wave number few nm. In 

this paper we employ spin-singlet isotropic s-wave pair¬ 
ing for simplicity to grasp the essential features of the 
PPE. Among the orbital and spin parts of the pairing 
symmetry, essential assumption in the present theory is 
that the PPE works in the spin part of the pairing. In 
addition to the case of the spin-singlet pairing, we can 
expect similar behavior of PPE also in the spin-triplet 
pairing case if the d-vector has components in the ah 
plane. The assumption for the orbital part as isotropic 
5-wave is not intrinsic condition. Both the 5-wave and 
the d-wave pairing show similar high field behaviors of 
PPE [45]. Also in the chiral p-wave pairing, we see the 
similar transverse components of the internal fields [46] . 
Thus, the replacement of the orbital part from the 5-wave 
pairing to the chiral p-wave pairing is possible, and we 
expect similar behaviors there, if the PPE is active in the 
spin-part of the pairing. 

The arrangement of this paper is as follows. In Sec. [n| 
we give the formulation based on Eilenberger theory with 
PPE. The spatial structures of vortices, including inter¬ 
nal magnetic field B(r) and the paramagnetic moment 

The form factors re- 
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^para(i*), are described in Sec. 
sponsible for SANS experiments are evaluated both for 
the longitudinal and transverse components relative to 


the applied field orientation in Sec. IV In the next sec¬ 
tion [V| we calculate the magnetization curves and Som- 
merfeld coefficient 7(5) as a function of magnetic fields 
to examine the first order transition’s jumps of these 
quantities at idc 2 - The distributions of P{B) of B{t) 
and P{M) of Mpara(r) of the vortex lattice state, that 
are responsible for the resonance line shape of the NMR 


spectra, are calculated in Sec. The magnetic torque 
curves are also evaluated in Sec. Eul We discuss in¬ 


trinsic anisotropy of Tvl and T ^^2 ia Sec. VIII The 


final section IX is devoted to conclusion and future prob¬ 
lems. The present paper belongs to our series of papers 
on the magnetic field orientation dependence of uniaxial 
superconductors: chiral p-wave case [46] and s-wave and 
d-wave cases without PPE HZ). 


II. QUASICLASSICAL THEORY INCLUDING 
PAULI PARAMAGNETIC EFFECT 


First, we explain the coordinate and the Fermi sur¬ 
face used in our calculations. We consider the case when 
the magnetic field orientation is tilted by 0 from the c 


axis towards the ah plane. We write the crystal co¬ 
ordinate as (a^h^c). To describe the vortex structure, 
we use the coordinate r = {x^y^z) where 2 : axis is set 
to the vortex line direction. Thus, the relation to the 
vortex coordinate and the crystal coordinate is given by 
{x^y^z) = {a^h cosd + csind, ccos 0 — hsinO). 

As a model of the Fermi surface, we use quasi-two 
dimensional Fermi surface with rippled cylinder shape. 
In the crystal coordinate, the Fermi velocity is assumed 
to be V = {va^Vb^Vc) oc (cos 0, sin0, sinpc) at p = 
{PaiPb^Pc) oc (pf cos 0 ,pf sin0,pc) on the Fermi sur¬ 
face [48|. From the Fermi surface, anisotropy ratio of 
the coherence lengths is estimated as 

r = ^c/Cb ~ (vc)p^V(^'b>y^ (1) 

where (• • •)p indicates an average over the Fermi surface. 

The spatial structure of quasiparticles in the super¬ 
conducting state is studied by the Eilenberger theory. 
Quasiclassical Green’s functions /(cc;^,p,r), /’^‘(co’n, p, r), 
and p(co’n, p, r) are calculated in the vortex lattice states 
by solving Riccati equation, which is derived from the 
Eilenberger equation 

{w„ + iM-B(r) + V ■ (V + iA(r))} / = A{r)g, 

K + - V ■ (V - iA(r))} p = A*{r)g, (2) 

in the clean limit with 

v-V 5 = A*(r)/-A(r)/t, (3) 

p = (1 — Rep > 0, and Matsubara frequency 

oJn im l48ll5Q] . The paramagnetic parameter p = 
Pb^o/tt/cbTc is proportional to the Maki parameter. We 
calculate the spatial structure of p in a fully self-consist 
way without using Pesch’s approximation m- We con¬ 
sider the case of isotropic 5-wave pairing, because the 
paramagnetic effect does not seriously depend on the 
pairing function of the orbital part [45]. Normalized 
Fermi velocity is v = v/uf with uf = We have 

scaled length, temperature, magnetic field, and energies 
in units of Tc, 5o, and tt/cbTc, respectively, where = 
hv^l^irkBPc and Bq = 0o/27r^o- is the flux quantum. 
The vector potential A = |B x r + a(r) is related to the 
internal field as B(r) = V x A = {Bx{r), By{r)^ Bz{r)) 
with B = (0,0,5), Bz{r) = B^hz{r) and {B^^By^hz) = 
V X a. The spatial averages of Bx^ By, and hz are zero. 
5 is the averaged flux density of the internal field. 

The pairing potential A(r) in the isotropic 5-wave pair¬ 
ing is calculated by the gap equation 

A{r)=goNoT ^ + (4) 

O^UJri 7^cut 

where po is the pairing interaction in the low-energy band 
\^n\ < and No is the density of states (DOS) at the 
Fermi energy in the normal state, po is defined by the 
cutoff energy Wc as (floAo)”^ = InT + 2TYZPo ojp. 
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We carry out calculations using the cutoff ujc = 
Current equation to obtain a is given by 

20 /cbTc. 


y/c 

V X V X a(r) = js(r) + V x Mpara(r) 

(5) , , 

(a) 

\ ^ 

\ ' ^ 

• 


with the screening current 

Y1 , (6) 

O^UJn, 

and the paramagnetic moment 

Mp,ra(r) = Mo ^ E (I“{5})pj -(7) 

Here, the normal state paramagnetic moment Mq = 
and hz = Hq/tt/cbTcV^StEVo. The Ginzburg- 
Landau (GL) parameter hz is the ratio of the penetration 
depth to coherence length for B || c. 

We set unit vectors of the vortex lattice as 

Ui = c(a/2, -\/3/2), U 2 = c{a/2, V3/2) (8) 

with = 2(j)o/{\/3aB) and a = 3T{d) |3S], as shown in 
Fig. f^a). We use the anisotropic ratio T{0) = ^y/^x ^ 

{vl)p /{vl)p^, that is, 


(b)|A(r)| 


(c) A^para(r) 


(d) 



(e) B,{r) 


m By{T) 



m = . , = (9) 

V cos^ 0 + F”^ sin^ 0 

of the effective mass model. Supposing the case of 
Sr 2 Ru 04 [1], we set to be k = 2.7 and the anisotropy 
ratio r(^ = 90°) = F = 60, which is suggested by the 
SANS experiment [4]. By the iteration of calculations by 
Eqs. 0 at T = O.lTc, we obtain self-consistent solu¬ 
tions of A^), A(r), and quasiclassical Green’s functions. 

From the selfconsistent solutions, we calculate the fol¬ 
lowing physical quantities. In Eilenberger theory, free 
energy is given by 


F = K2(|B(r)-B|2)r-M^(|B(r)|2)r 



(A/t + A*/) 



when Eqs. © and @ are satisfied [48]. (• • •)r indicates 
the spatial average. The magnetization is calculated as 
A^totai = B — where the external field H is given by 


FIG. 1: (Color online) (a) Unit vectors ui and U 2 of the 
vortex lattice. Circles indicate vortex centers. Gray region is 
a unit cell of our calculations, (b) |A(r)|. (c) Mpara(i*). (d) 
^^(r) — B. (e) Bx{r). (f) By{r). In (b)-(f), we show density 
plot within a unit cell, when 0 = 89° at H = 1.5 and /x = 0.04. 


Aftotai “ Afpara- As the resonaucc line shape of the NMR 
spectrum for the Knight shift, we calculate the distri¬ 
bution function P{M) = {6{M — Mpara(i*)))r from the 
spatial structure of Mpara(i*)- On the other hand, in 
the case of negligible hyperfine coupling, the NMR signal 
shows “Redfield pattern” given by the distribution func¬ 
tion P{B) = (d(H — H(r)))r calculated from the internal 
field B{r). 

When we calculate the electronic states, we solve Eq. 
© with iuOn Epir]. The local density of states (LDOS) 
is given by A(r, E) = A-^(r, E) + E), where 






k^B 


/ D / \T r 1 I / (/^^ + 1 

P { 5 } + ^-TITT-f 


( 

0<Cc>, 

PuJnRe{g — l})k) 


( 11 ) 


from Doria-Gubernatis-Rainer scaling [52l[53]. The para¬ 
magnetic and diamagnetic components of the magnetiza¬ 
tion are, respectively, Mpara = (Afpara(i*))r and Mdia = 


N^{r,E) = No{Re{g{uJn P 

with a = 1 (—1) for up (down) spin component. We 
typically use g = 0.01. The DOS is obtained by the spa¬ 
tial average of the LDOS as N{E) = N^{E) P A^(E) = 
(A(r, £'))r. We consider the H-dependence of the Som- 
merfeld coefficient of the specific heat given by the zero- 
energy DOS j{B) = N{E = 0)/Ao, and the paramag¬ 
netic susceptibility Xspin(5) = (Mpara(r))r/Afo. These 
are normalized by the normal state values. 


















4 



FIG. 2: (Color online) (a) dependence of the pair potential 
when 0 = 89° (solid lines) and 0 = 90° (dashed lines). Spatial 
averaged values of |A(r)| are presented for /x = 0 and fi = 
0.04. The latter (/x = 0.04) exhibits first order transitions for 
both 89° and 90°. 


III. SPATIAL STRUCTURES OF VORTICES 


To discuss 5-dependence of the internal field distribu¬ 
tion B(r) = V X A, we consider flux line lattice (FLL) 
form factors F{cih,k) = iFa:(h,k), Fy(h,k), F^{h,k)), which 
is obtained by Fourier transformation of the internal 
field distribution as B(r) = F(q/i,fc) exp(iq/i,fc ■ r) 

with wave vector = hqi + kq 2 - h and k are in¬ 

tegers. Unit vectors in reciprocal space are given by 
qi = { 2 TT/c){l/a,-l/V 3 ) and q2 = { 2 tt/ c){l/a,l/\/Z). 
The 2:-component from Bz(r) gives the inten- 

sity of conventional non-spin flip SANS. The transverse 
component, |Ftr(h,k)l^ = \Fx(h,k)? + is accessi- 

ble by spin-flip SANS experiments [4l 154]. 

First, we study the vortex states when the magnetic 
field orientation is tilted by 1° from the ah plane {0 = 
89°). In Fig. we show the calculated spatial struc¬ 
tures within a unit cell of vortex lattice at B = 1.5 and 
/i = 0.04. The pair potential A(r) has phase winding 
27r at the vortex center, and the amplitude |A(r)| in Fig. 
[^b) is suppressed at the vortex core. The paramagnetic 
moment Mpara(i*) Fig. j^c) is suppressed outside of 
the vortex core. Mpara(r) appears within the narrow re¬ 
gion around the vortex core and has a large peak at the 
vortex center. As shown in Fig. Eld), the ^-component 
of internal field, has a peak at the vortex cen¬ 

ter, and decreases as a function of the distance from the 
center. The peak height of Bz{r) is enhanced by the 
contribution of Mpara(i*) at the vortex core [45]. The 
vortex state has a conventional spatial structure in the 
vortex lattice also when B is tilted from the ab plane, 
if length is scaled by the effective coherence length in 
each direction. When B is tilted from the ab plane, the 
transverse components Bx{t) and 5^(r) appear in the in¬ 
ternal field distribution, as shown in Figs. [^e) and Bf). 
The magnitude of 5^(r) is larger than that of Bx{r). 
The stream lines of By{r) in Fig. 00 flow towards —y 
direction along vertical stripe region connecting vortex 
cores. Between the neighbor stripe regions, the stream 



FIG. 3: (Golor online) 5-dependence of the FLL form fac¬ 
tors. (a) |52 (i,o)|^ and |52(i,i)|^ for /x = 0, when 0 = 89° 
(solid lines) and 90° (dashed lines), (b) The same as (a), but 
for fi — 0.04. (c) |5a,(i 0)1^ for /x = 0 and 0.04 when 6 — 89°. 
(d) |5^(i 1)1^ for /X = 0 and 0.04 when 0 — 89°. 


line flows towards +^-direction. The weak contribution 
of 5a, (r) in Fig. [^e) indicates that the stream lines have 
weak counter-clock-wise (clock-wise) winding at positive- 
X (negative-x) region near vortex core. These stream line 
structures of the transverse field is qualitatively the same 
as those obtained by London theory m, and as those in 
a chiral p-wave pairing [46] . 

When the paramagnetic effect is not considered (/x = 
0), the upper-critical field is Hc 2 ,c = 0.56 for B || c and 
Hc 2 ,ab = 43 for B || ab, reflecting large anisotropy F. 
Figure [^presents the amplitude of the pair potential as 
a function of 5. In the case /x = 0.04, the paramagnetic 
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pair-breaking is negligible for B || c so that i^c2,c is un¬ 
changed. However, for B || a6, the paramagnetic pair 
breaking becomes eminent at high fields and limits the 
upper critical field to Hc 2 ,ab = 9.1. The phase transi¬ 
tion at Hc2,ab becomes first order as coinciding with the 
observation in Sr 2 Ru 04 at low temperatures [HE]. 

For the field orientation tilted by 1° away from the ah 
plane, namely 0 = 89°, Hc2 is suppressed from Hc2,ab 
at 0 = 90° as seen from Fig. It is noted that those 
Hc 2 suppressions are quite different: While in the /i = 
0 case Hc2{0 = 89°)/i7c2(<9 = 90°) = 32/43 - 0.74, 
approximately satisfying the expectation based on our 
effective mass model; r{0 = 89°)/r(6> = 90°) ^ 0.69, the 
Hc 2 suppression in the /i = 0.04 case is very small and 
remains first order. This is because Hc 2 is determined 
by the PPE and controlled by the Pauli paramagnetic 
critical field Hp{0) which has a weak 0 dependence [45] . 
This point will be discussed later in connection with the 
nature of the phase transitions. 


IV. FLUX LINE LATTICE FORM FACTORS 
A. Longitudinal component 

We discuss the 5-dependence of the FLL form factor 
for 9 = 90° and 89°. Figures |^a) and ^b) present the 
5-dependence of |T 2 (ip)P and |T^(i,o) • These corre¬ 
spond to the intensity of the non-spin-flip SANS experi¬ 
ments. When /i = 0 in Fig. |^a), |T^(i,i)P and |52(i,o)P 
show exponential dec^ as a function of 5. However, 
when /i = 0.04 in Fig. j^b), both |T^(i,o)P and |5z(i,i)P 
increase towards a maximum at 5 ^ 7 for 6> = 90°. 
This increasing behavior is due to the enhancement of 
the paramagnetic moment at the vortex, which enhances 
the peak of Bz{r). This mechanism [45] was discussed 
to explain the 5-dependence of the SANS intensity in 
CeCoIns [12] [56] and TmNi 2 B 2 C [57] • Compared to the 
case of 0 = 90°, the intensities of \Fz(h,k)\‘^ are more en¬ 
hanced for 0 = 89° at low fields. This is because the 
intensity \Fz(h,k)\‘^ is roughly related to the effective GL 
parameter kq as \Fz(^h,k)\‘^ oc By the anisotropy of 

V in Eq. K,o ^ bvr{0) for the field orientation 0. 

Thus, KgQo < KgQo. At high fields, the peak position of 
\Fz{h,k)\‘^ is shifted to 5 ^ 6 when 0 = 89°, reflecting the 
6>-dependence of anisotropic i7c2- 


B. Transverse components 

The 5-dependence of the transverse component is 
shown in Eigs. |^c) and |^d). As for (1,0) spot, 
l^’.d.oTT Since |F,(i,o)|" < lO/ii. 
decreases monotonically as a function of B. As for (1,1) 
spot, |Ftr(i,i)P ~ since -- 0. |i^s/(i,i)f 

decreases as a function of H, after it increases at low B. 
As in the chiral p-wave pairing [46], |F^(i,i)P has large 



FIG. 4: (Color online) 5-dependence of the FLL form factors 
when 0 — 89°. We plot renormalized values |F^(i i)/i7c2,a6|^ 
as a function of for \i — 0.04 and 0. The points -h 

indicate experimental values [3] on Sr 2 Ru 04 . We also present 
\Fz{i,q) / Hc2,ah\^ and \Fz{i^i)/ Hc2,ah\^ ■ The vertical axis is a 
logarithmic scale. 


intensity, compared with |5c(i^o)P \Fz{h,k)\‘^- This is 

consistent to the fact that only the spin-flip scattering at 
(1,1) is observed in the SANS experiment [4] on Sr 2 Ru 04 . 
From Figs. [^c) and [^d), we see that the enhancement 
due to the paramagnetic effect does not appear in the 
transverse component p. Rather de¬ 

creases rapidly at higher fields, reflecting the paramag¬ 
netic suppression of superconductivity. 

For the quantitative comparison with the experimental 
data [4] in Sr 2 Ru 04 , we discuss the form factors and 5 
in unit of Hc 2 ,ab as plotted in Fig. In the case /i = 
0.04, i)/i7c2,a6p is larger because Hc2,ab is smaller. 

In Fig. we also show the experimental data [4] on 
Sr 2 Ru 04 with Hc 2 ,ab = 1-5[T], i.e., Ftr(i^i) = 0.677[mT] 
at 5 = 0.5[T] and Etqi,!) = 0.485[mT] at 5 = 0.7[T]. 
The magnitude of i)/i7c2,a6p in experimental data 
can be quantitatively reproduced in the case p = 0.04 
including the effect of B[c2,ab suppression. From Fig. 
we also see that 10“^-times finer resolution is necessary 
in the SANS experiment to observe the spot of \Fz(^h,k)\‘^ 
for the non-spin-flip scattering, which is expected to be 
an increasing function of 5 at the middle field range. 

The ^-dependence of the |F^(i,i)P is presented in Fig. 
1^ As a function of 6>, |F^(i,i)P increases until a peak 
near 90°. After the peak it decreases rapidly towards 
zero just at 90°. At low enough field 5 = 1.5, |F’^(i,i)P 
shows similar behavior both for /i = 0 and 0.04. With 
increasing 5, the peak position is shifted to higher 6>, 
and the amplitude is decreased. At higher fields such as 
5 = 13.5, |F^(i,i)P becomes very small for /i = 0. How¬ 
ever, these high field regions vanish for /i = 0.04 because 
of high-field suppression of superconductivity. For the 
quantitative comparison. Fig. [^c) shows renormalized 
values i)/i7c2,a6p in n logarithmic scale with the 
SANS results on Sr 2 Ru 04 for two cases B/H^ 2 ,ab ^ 0-33 
(5 = 3 .OH 0 and 0.5[T]) and B/H^ 2 ,ab ^ 0.5 {B = 4.5Ho 
and 0.7[T]). We see that the experimental data are well 
fit by the theory for /a = 0.04 near 0 = 89°. The theoret¬ 
ical values for /i = 0 is very small compared to the SANS 
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FIG. 5: (Color online) Field orientation ^-dependence of 
the transverse FLL form factor, (a) \Fy(^i^i^/Hc 2 ,ab\‘^ as a 
function of ^ for /x = 0 at B = 1.5, 4.5, 9.0, and 13.5. (b) 
i)/Hc 2 ,ab\‘^ as a function of ^ for /x = 0.04 at B = 1.5, 
3.0, and 4.5. (c) \Fy(^i^ i)/Hc 2 ,ab\^ in a logarithmic scale as a 
function of 0 for /x = 0.04 at B/Hc 2 ,ab ~ 0.33 and 0.5 {B = 3.0 
and 4.5). The points + indicate experimental values [4] on 
Sr 2 Ru 04 at 0.5[T] and 0.7[T]. We also plot \Fy(^i^i^/Hc 2 ,ab\‘^ 
for /X = 0 at R = 9.0 {BlHc 2 ,ab = 0.21). (d) 0 dependence of 
the longitudinal FLL form factor \Fz(i^i)/Hc 2 ,ab\‘^ for /x = 0.04 
at BlHc 2 ,ab = 0.17, 0.33, and 0.50. 


results. 

We also plot the longitudinal component 
\^z{i,i)/H c2, ab\‘^ for = 0-04 in Fig. j^d) that is not yet 
observed in Sr 2 Ru 04 . It is seen that |F^(i i)/iLc 2 ,a 6 p 
grows as B increases at ^ = 90° because of PPE, as 
seen in Fig. lib). At a low field B/H^ 2 ,ab— 0.17, 
|^z(i,i)P monotonically increases when 0 decreases, 
since 1^2(1,!)P oc We notice that the longitudinal 

components of the form factor are already observed 


for 0 = 0, i.e., H \\ c [58]. The detailed analysis of 
those form factors has not done yet, but it seems to 
be similar to the results for the square lattice for the 
d-wave case [59]. For larger fields B/Hc 2 ,ab= 0.33 and 
0.50, it takes a peak at finite 0 because the effective 
magnetic field 5 /iLc 2 becomes large as 0 decreases 
from 90°. Then ^^( 1 , 1 ) vanishes ultimately towards Hc 2 
where the order parameter is zero. This peak behavior 
in is similar to those shown in Fig. j^b) where 

the longitudinal components as a function of B exhibit 
peaks just below iLc 2 - 


V. JUMPS AT FIRST ORDER Hc2 TRANSITION 
A. Magnetization curves 

We calculate the magnetization curves for 0 = 89° and 
90° both in the cases of /x = 0 and /x = 0.04 at T = 
O.lTc as shown in Figs, [^a) andj^b). In the /x = 0 
case, Mtotai(5) corresponds to that of an ordinary type 
II superconductor, because Mtotai(5) comes exclusively 
from the orbital diamagnetism due to the orbital current. 
Since the second order transition occurs at Hc 2 in this 
case, Mtotai(5) smoothly becomes zero. 

As seen from Fig. |^b) in /x = 0.04, Mtota\{S) exhibits 
the jumps AMtotai at Hc 2 both for 0 = 89° and 90°, 
corresponding to the first order transition. The magne¬ 
tization jump AMtotai consists of the two components; 
the orbital diamagnetism AMdia and the spin paramag¬ 
netism AMpara- ^ = 90° Mpara = 1-3 X 10“^ and 
^dia = —0.7 X 10“^ at 5 = Hc 2 = 9.1. Thus AMtotai is 
dominated by the spin paramagnetic component. Since 
at 5 = iLc 2 , Mnormai = 2.0 X 10“^, the relative jump 
AMtotaiZ-M'normai = 38.5%. As Seen from Fig. |^b), the 
0 = 89° case also gives rise to a similar AMtotai value. 

This number is favorably compared with the exper¬ 
imental value AMtotai/^normal = 15% in Sr 2 Ru 04 

at low temperatures [8]. A slightly larger value of 
AMtotai/^normal in our Calculation can be remedied by 
considering the multiband effect because near Hc 2 the 
minor band may be almost in the normal state where 
the minor gap is already vanishing. The minor band 
contribution can be estimated as AMtotai/(ATnormai + 
ATminor) ^ 19% bccausc the DOS of the minor band is a 
half of the total DOS. 

As seen from Fig. [^c), the contribution of the or¬ 
bital diamagnetism Mdia to the first order jump amounts 
to AMdia = —0.7 X 10“^. The weight of the jump, 
AAfdia/Mnormai = 3.5%, is an ordcr too smaller than the 
observed value. Thus without PPE it is impossible to 
understand the large magnetization jump. We also point 
out that the magnetization curve for the chiral p-wave 
case (see Eig. [^a) in Ref. [46|) is almost same as in usual 
type II superconductor without PPE shown in Eig. [^a). 
Thus, if we assume a hypothetical first order transition 
at Mist, then Mist ^ 0.25Mc2 to account for the magneti¬ 
zation jump AMtotai/Afnormai = 15%. So far there is no 
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known theory to explain the first order transition in the 
chiral p-wave pairing. 


B. Specific heat jump at Hc2 

We show the calculated results of the DOS at the 
Fermi level at low temperature T = O.lTc in Fig. [^d), 
which corresponds to the Sommerfeld coefficient ^5), 
namely C/T at low temperatures in the superconduct¬ 
ing state. It is known by the explicit calculations [45] 
that 7(5) is approximately scaled to the spin suscepti¬ 
bility Xspin(5) as is seen from Fig. [^d). This is because 
both quantities Xspin(5) and 7(5) come from the same 
DOS of the quasi-particles near the Fermi level [45|. In 
fact the experimental value [7] of the specific heat jump 
at Hc 2 is Aq/qnormai ^ 10%, roughly coinciding with 
AMtotai/^normal = 15% mentioned above. As is seen 
from Fig. |^d), the jump of "f{B) is slightly smaller than 
that of Xspin because two quantities are not exactly iden¬ 
tical where the former is an integration of the DOS over 
/j^bB while the latter is the DOS at the Fermi level. 

It should be noted that the entropy jump [6] probed by 
the magneto-caloric measurement is consistent with the 
specific heat jump as discussed in Ref. [7], meaning that 
three experiments, magneto-caloric, specific heat and 
magnetization are mutually consistent with each other. 
If this identification is true, the Knight shift should de¬ 
crease as shown in Fig. |^d), which is contrasted with the 
claim by the NMR experiments [TSHTS] where the Knight 
shift remains unchanged, irrespective to nuclear species 
(i7o, ^^^Ru, and ^^Ru), the field orientations and 

the field values. This is quite at odds in the present anal¬ 
ysis. 


VI. FIELD DISTRIBUTIONS 

Figures [^a) and [^b) display the field evolutions of 
P{B) and P{M) together with the contour maps of Bz{r) 
and Mpara(i*) within a unit cell. It is seen that as in¬ 
creasing field towards Hc 2 the vortex core site and its 
surrounding sites exclusively accommodate the param¬ 
agnetic moments induced by PPE where the highest Bz 
and Mpara are situated. The mean value of P{M) equals 
Xspin in Fig. I^d). At 5=8.5, due to the contributions of 
the paramagnetic moment enhanced at the vortex core, 
P{B) and P{M) have larger weights near the small peak 
at highest edges, whose positions of highest edge exceed 
5 and Mo, respectively. Thus the so-called Redfield pat¬ 
tern P{B) is strongly modified from the standard asym¬ 
metric distribution in ordinary superconductors (60] |6T] , 
such as Nb [62]. This is also true for P(M) where the 
asymmetric pattern is modified so that the higher M 
range of the spectrum is enhanced. 

Those asymmetric spectrum patterns should be ob¬ 
served by NMR experiments, where neither asymmetric 
P{B) nor P{M) patterns are not observed in any nuclear 



FIG. 6: (Color online) 5-dependence of the magnetization, 
(a) Mtotai(5) = Mdia(5) for /x = 0, when 0 = 89° (solid 
line) and 90° (dashed line). In the inset Mtotai(5) is enlarged 
near 5c2. (b) Mtotai(5), Mpara(5) and Mdia(5) for fi = 0.04 
when 0 — 89° (blue solid lines) and 90° (red dashed lines), 
(c) Mdia(5) is focused near Hc 2 for 0 = 89° and 90° to see 
the jumps of AMdia at 5c2. (d) The scaling behaviors of 

7(5) and Xspin(5) for 0 = 90°. The jumps of Ay and Ayspin 
relative to its normal values are seen at the first order 5c2 
transition. 


species (^^O, ^^Sr, ^^^Ru, and ^^Ru ) for B || ah. They re¬ 
main the same patterns as in the normal state [T3|. Note 
that the characteristics in Fig. are indeed observed 
in CeCoIns m- On the other hand, for B || c, a clear 
Redfield pattern is observed by the muon spin resonance 
experiment |63| . By analyzing this pattern they correctly 
reduce the vortex lattice symmetry, namely a square lat¬ 
tice that is confirmed later by SANS experiments [58] . 
















FIG. 7: (Color online) Topographic maps of (a) ^^(r) — B 
and (b) Mpara(r) within one unit cell at B = 1.5, 4.5 and 8.5 
for fi = 0.04 and 0 = 90°. The held distribution (c) P{B) and 
(d) P{M) associated to (a) and (b), respectively. 


VII. MAGNETIC TORQUE 

Since we obtained the self-consistent solutions of Eilen- 
berger equation under a given T and 5, it is not difficult 
to calculate the magnetic torque t{0) = dF/dO by us¬ 
ing the free energy F as a function of 0. The obtained 
free energy F{0) is displayed in Figs, [^a) andj^b) for 
/i = 0 and fi = 0.04, respectively. It is seen that for both 
cases all the free energy curves smoothly become zero 
when 0 decreases away from 0 = 90°, meaning that those 
are all second order Hc 2 transitions in the held range 
5<Fc2(<9 = 88°). 

Figures l^c) andj^d) show the magnetic torque curves 
t{0) for /i = 0 and /i = 0.04, respectively. It is seen from 
those that the sharp minima in r{0) for both cases are 
located just near 0 = 90°. The fact that the minimum 
position Omin in is conhned near 0 = 90° is due to 
the large uniaxial anisotropy F = 60. This behavior is 
easily ht by the Kogan torque formula [64] based on the 
London theory: 

sin2<9 f]THc2,c 

r{0) oc = In ^ ^ 

V cos^ 0 sin^ 0 By cos^ ^ -f F^ sin^ 0 

(13) 

with 0 = 90° — where r) is a coefficient with the or¬ 
der ^ 1. The minimum ^min occurs at ^min = 88.7° for 
F = 60 with T] = 1.5, which is consistent with our Eilen- 
berger solution. It should be noticed that at lower fields 
t{0) is insensitive of the presence or absence of PPE ac¬ 
cording to our results in Eigs. ic) andj^d). Thus, both 


cases are described by the Kogan formula which only 
depends on F. In fact the minima observed experimen¬ 
tally show r\j 89° at higher fields, which becomes 
towards lower fields (see Fig. 3(d) by Ref. 
Ej). Also in the numerical calculation with F = 60, ^min 
shows similar weak 5-dependence. Thus the torque data 
support the large uniaxial anisotropy with F = 60 for 
Sr 2 Ru 04 . We note that if we choose F = 20 as indicated 
by Hc2,ab/Hc2,c r\j 20, we find ^min ^ 87° which is far off 
the experimental datajS] within the present experimen¬ 
tal accuracy. Thus the intrinsic anisotropy of Sr 2 Ru 04 
should be F = 60 rather than F = 20. The latter number 
is now understood as arising from the suppressed Hc 2 by 
PPE. 

Since the magnetic torque is related to the transverse 
components of the internal field, it is interesting to com¬ 
pare the \t{ 0 )\ curves with the form factor |F^i for 
both /i = 0 and ji = 0.04 as shown in Eigs. ^ (a) and 
ib). An approximate scaling relationship between them 
is seen from both cases. In particular, the maximum po¬ 
sition ^max in both quantities yields the same value for 
the higher field data. This is indeed seen experimentally 
(see Eigs. 3(c) and 3(d) in Ref. |8|). 

In Eig. [Tol^a) we compare the theoretical torque curves 
and corresponding experimental data [8] for selected field 
values. It is seen from it that they show a good agree¬ 
ment, in particular in the higher field data, including 
the maximum angles and vanishing angles of the torque 
curves. The highest field theoretical curve 5=8.6 ex¬ 
hibits a first order jump at 0 = 89.1°, which nicely co¬ 
incides with the data at 1.4[T]. On the other hand, the 
lower field data at 0.2 [T] show a deviation from the theo¬ 
retical curve 5=1.5 in their maximum angles. This may 
come from the multi band effect, which will be discussed 
in the forthcoming paper [65] . 

In Eig. [Toj^b) we summarize our maximum angle 
data coming from the torque curves and the transverse 
form factors and compare those with the experimental 
data m- As already indicated in EigJ^ the discrepan¬ 
cies of the maximum angles between the torque and form 
factor occur when the field is lowered. Since the form fac¬ 
tor data at lower fields is lacking at present, we cannot 
judge whether or not those discrepancies are strength¬ 
ened further by future SANS experiments. Except for 
those lowest field data the overall agreement seems to 
be satisfactory. In other words, the present single band 
theory gives a reasonable explanation to those data. 


VIII. DISCUSSIONS ON PHASE DIAGRAM 
AND INTRINSIC ANISOTROPY 

In previous sections, the vortex lattice anisotropy 
rvL(^) = <^/3 in the definition of Eq. § is assumed to 
be given by the effective mass model in Eq. <§• We also 
perform calculations to determine Fvl(^) by the min¬ 
imization procedure of the free energy, which is much 
time consuming process compared above. The results 
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FIG. 8: (Color online) The ^-dependences of the free ener¬ 
gies for various helds B. (a) /i = 0 and (b) fi = 0.04. The 
corresponding torque curves t{0) — dF/dO. (c) /x = 0 and (d) 
/X = 0.04. 


are shown in Fig. pT[a). It is seen that this yields a 
slightly larger Fvl (^compared with the effective mass 
model shown by a line there around 0 = 90° ± 2° region, 
beyond which all data points tend to coincide with a line 
of the effective mass model. We confirm that this de¬ 
viation of Fvl(^) does not alter our results in previous 
sections in a serious way. 

We note that, as presented in Fig. pTj^a), experimen¬ 
tal data [4] also slightly deviate from the effective mass 
model for 87° < 0 < 89°, which is similar to Fvl(^) by 
the free energy minimum. This behavior will be discussed 
in a forth coming paper based on multiband model [65] . 
From the ^-dependence of Fvl(^) in Fig- [Tit a), Ibe in- 



FIG. 9: (Color online) The scaling behaviors between 

|Fy(l,l) I and \r{0)\ as a function of 6. (a) /x = 0 and (b) 
= 0.04. Magnitude of each quantity is scaled by the maxi¬ 
mum value. 


trinsic uniaxial anisotropy of the system can be identified 
as F = 60. This number just corresponds to the Fermi ve¬ 
locity anisotropy of the p band, since the band-dependent 
anisotropies are estimated as Ta = 117, F /3 = 57, and 
F^ = 174 for the a, /3, and 7 Fermi-surface sheets, re¬ 
spectively, according to the dHvA experiments [T]. We 
emphasize that this is not accidental, if the [3 band plays 
a major role to govern the mixed state properties of the 
total system in high fields, further suggesting that the 7 
band plays a secondary role, which is contrary to what 
many previous works claim, such as Nomura and Ya- 
mada [ 20 ] . 

According to the present analysis, the Hc 2 ,ab suppres¬ 
sion is explained by PPE. We evaluate Hc2{0) for each 
0 by the estimate of the critical field where the order 
parameters vanish on raising 5, as done in Fig. 

The ^-dependence of Hc2{0)/Hc2c is presented in Figs. 
[Fib) and m c), where the filled (empty) symbols corre¬ 
spond to the first (second) order transitions. Our calcula¬ 
tion shows that the first order transitions only occur for 
0 = 89°, 90°, 91°, beyond which all Hc 2 transitions be¬ 
come second order as displayed in Figs. [TT]^b) and[IT]^c). 
This is consistent to experimental data^ presented in 


the figures. There, in Sr 2 Ru 04 the first order transitions 
occur within 90° ± 2°. The first order transition near 
the ab plane appears because th e eff ective paramagnetic 
parameter fl{0) = /iF( 6 >) in Eq. (A 2 ) exceeds the critical 
value /Xcr ^ 1-7 for the first order transition only for the 
angles 89° < 0 < 91°. 

[h) and 


As seen from Eigs. 
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c), Hc2{0)/Hc2,c val¬ 


ues by the numerical calculations are well fit by a simple 


function [ 66 ] in Eq. ( |A3| ) with /x = 0.04, explained in 
Appendix. The values of Hc2{0)/Hc2,c are slightly un- 
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FIG. 10: (Color online) (a) Comparison with the theoreti¬ 
cal torque curves \t{0)\ for 5 = 1.5, 3.0, 4.5, 8.6 (lines) and 
experimental data [8] for B = 0.2, 0.5, 0.7, 1.4[T] (empty sym¬ 
bols). We have adjusted the maximum values of the torque 
curves and displayed those curves by changing its maximum 
values arbitrarily to be clearly seen. The highest field theoret¬ 
ical data B = 8.6 and experimental data 1.4[T] clearly show 
the jumps associated with the first order transition, (b) Max¬ 
imum angles of the form factors (triangles) and the torque 
curves (circles) in the B and angle 0 plane. Theoretical re¬ 
sults (filled symbols) are compared with the corresponding 
experimental data [U [8] (empty symbols). In the scale of ver¬ 
tical axis, Hc 2 ,ab = 9-1 in theoretical estimate is assigned to 
be 1.5[T]. 




der the experimental values. This can be easily remedied 
by changing the /r value. Namely, instead of the present 
value /i = 0.04, the refined value (i = 0.0293 sho ws m uch 
better fitting to the experimental data by Eq. (A3), as 
shown in Figs. [^b) and[^c). 

We also evaluate the temperature dependence of the 
ratio Th^^{T) = i^c 2 ,ab(^)/^c 2 ,c(^), and compare it 
with the experimental data [67] in Fig. [TT]^d). Near T = 
Tc, both in numerical and experimental data, T 
shows the large anisotropy ratio, tending to ^ 60, which 
is governed by the Fermi velocity anisotropy ratio of the 
(3 band F^ = 57. Upon decreasing T, this ratio progres¬ 
sively becomes small because of the PPE. This is cap¬ 
tured by our numerical calculation, and the tending limit 
towards the low est T is 16. Th e T-dependence is well fit¬ 
ted by Eq. ( |A3| ) with Eq. (A4) for fi = 0.04, as shown in 
Fig. [TT](d). In the experimental data which reduces 20 


FIG. 11: (Color online) (a) The ^-dependences of the vortex 
lattice anisotropy Fvl(^). Open circles indicate the experi¬ 
mental data for B=0.5[T] [4]. Other symbols are for Fvl(^) 
evaluated by the free energy minimum at B jB q =2, 4, 8, and 
16 for fi = 0. The line presents F(^) of the effective mass 
model in Eq. § with F = 60. (b) The ^-dependences of 
Hc2{0)!Hc 2 ,c at T = O.lTc for /i — 0.04. (c) Enlarged figure of 
(b). The results of numerical calculations by Eilenberger the¬ 
ory are presented by square points. The experimental data [8] 
are shown by circles. There, the filled (empty) symbols indi¬ 
cate the first (second) order transition. The solid line shows 
F(^) of the effective mass model with F = 60. The dashed 
lines correspond to the theoretical curves calculated by Eq. 
( |A3| with F = 60, where ji — 0.04 for fitting to numeri¬ 
cal calculations and ji = 0.0293 for fitting to experimental 
data, (d) The anisotropy Fif^ 2 (^) = Hc 2 ,ah{T)/Hc 2 ,c{T) as a 
function of T. The experimental data[67] (circles) and the nu¬ 
merical results by Eilenberger theory with fi — 0.04 (squares) 
are shown. The three continuous lines are evaluated by Eqs. 
( A3t and (j^ with (1) F = 180 and /i=0.0293, (2) F = 60 
and ii — 0.0293, and (3) F = 60 and fi — 0.04. 
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at low T [67] is fitted by Eq. (A3) with /i=0.0293 as in 
the case of Figs. inib) and [TT](c). We also note that the 
fitting line with E = 180 largely deviates from the exper¬ 
imental data. A similar analysis on the Hc 2 anisotropy 
data m is performed by Choi |36| to come to the same 
conclusion. We point out again that in numerical calcula¬ 
tion (blue square) the point in low T range in Fig. [TT]^d) 
corresponds to the first order Hc 2 ,ab transition, while at 
least above T/Tc > 0.5, Hc 2 ,ab is of second order. Ac¬ 
cording to the experiments mil], the first order line at 
Hc 2 ,ah{T) extends to around T/Tc > 0.4 ^ 0.5. The 
accurate termination point between the first and second 
order transitions will be a future problem. 

In summary of this section, in both estimations 
of rvL(^) and the intrinsic anisotropy F of 

Sr 2 Ru 04 is identified as F = 60. Since it corresponds 
to = 57 given by the P band among the known three 
bands. The anisotropy F = 60 indicates that the P band 
is fully responsible for determining various observables. 
Thus this should be the major band, while the 7 band 
with = 174 is not appropriate for the major band 
as seen from Fig. [TT]^d) and must be the minority band 
and the a band plays a negligible role because its DOS is 
10% of the total. Those considerations partly justify the 
present single band model to grasp the essential points. 


spin-triplet pairing. Even in the spin-triplet pairing, the 
PPE works if the d-vector has a component within the 
ab plane. The essential characteristics of the mixed state 
magnetic properties investigated here remains intact even 
for those pairings. In the present paper we have taken 
the isotropic 5 -wave just for computational convenience 
and for the illustrative purpose. 

From the results that the PPE plays significant contri¬ 
butions on the superconducting properties for B || a6, we 
recognize the importance of the further studies to deter¬ 
mine the spin part of the pairing function. In Sr 2 Ru 04 , 
due to the spin-orbit coupling of the electronic states, we 
have to consider the coupling of the spin and the orbital 
parts in the spin-triplet pairing function [2^[68H7Q] . This 
may be an essential factor to realize a spin-triplet super¬ 
conductor where the d-vector has the component in the 
ab plane. 

Several important experimental facts remain unan¬ 
swered in the present theory: 

(1) As shown in Fig. [^c), the transverse form factors 
as a function of 0 under a fixed field in the SANS exper¬ 
iment |4] exhibit a strong decay of its amplitudes long 
before reaching the expected idc 2 , when 0 decreases away 
from the ab plane at d = 90° in the lower fields. We an¬ 
ticipate the interplay of the multiband effect to explain 


IX. CONCLUSION AND UNSOLVED 
PROBLEMS 

The essential assumption in the present theory is that 
the PPE works in the spin part of the pairing function. 
The assumption for the orbital part as isotropic s-wave is 
not intrinsic in our calculations. As the pairing function 
inducing the PPE, in addition to the spin-singlet pairing, 
the spin-triplet pairing is also available if the d-vector has 
components in the ab plane. 

There exist several outstanding experiments to claim 
as the evidence for a spin-triplet chiral p-wave pairing re¬ 
alizing in Sr 2 Ru 04 . Among them, only the Knight shift 
experiments USHIS] by NMR and polarized neutron scat¬ 
tering are treated as an evidence for that the spin part of 
the pairing function is the spin-triplet. There the Knight 
shift remains unchanged for any magnetic field orienta¬ 
tions, any id, T and any nuclear species available so far. 
These are against the present theory. However, the un¬ 
changed T-dependence for all field orientations is a mys¬ 
tery, since the Knight shift is expected to be decreased for 
some of field orientations even in the spin-triplet pairing. 

Other experiments than the Knight shift are evidences 
for the orbital part of the pairing. There, the time re¬ 
versal symmetry breaking is suggested by muon spin ro¬ 
tation (/iSR) |40|, which shows the appearance of the 
spontaneous moments below T^. These varieties of the 
orbital part can be easily accommodated in the present 
theoretical framework. As for the time reversal sym¬ 
metry breaking, the possible pairing can be an 5 + id 
or d + id in the spin-singlet pairing, or px + ipy in the 


( 2 ) It is also obvious that the 7 (id) behavior shown in 
Fig. [^d) is quite different from the observed C/T at 
low T especially at lower field regions m where a ^/H- 
like sharp increase is observed. This was interpreted as 
coming from the minor band (a + /3) that is assumed to 
have a line node like gap structure simply because the 
DOS value attained in that field region seems to that for 
the a-\-P bands (43% of the total DOS). The 7 band with 
57% of the total DOS is regarded as the major band. 

This interpretation is at odds with the present the¬ 
ory in several points; If this is true, the “major” 7 band 
should be responsible for the high field region. Since the 
Fermi velocity anisotropy Tj = 174 for the 7 band, we 
would expect Fvl ^ 174 for B || ab. There is no indica¬ 
tion for it in the available SANS data |4] and other EZ]. 
Instead, the SANS data a tvl r\j 60 indicate that the 
high field phase should be the /3 band. We point out also 
the data [67] iFc 2 ,ab(^)/^c 2 ,c(^) shown in Fig. lin'i) 
which is directly related to the Fermi velocity anisotropy 
in Eq. 0 near Tc. Thus the P band is also respon¬ 
sible for it [65]. There is no trace in the existing data 
to show that the the 7 band plays a major role. The 
present single band theory assumes the (3 band as the 
major band, neglecting the minor 7 band. We need to 
refine it by taking into account both bands in addition 
to the remaining a band. A multiband theory based on 
Eilenberger framework belongs to a future work [65] . 

In order to further advance the Sr 2 Ru 04 problem con¬ 
cerning its pairing symmetry and multiband nature, we 
propose here several decisive experiments: 

(A) The SANS experiments to observe the longitudinal 
component Fz in the FLL form factors. As already pre- 
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dieted in Figs. and |^c) the magnitudes i^z(i,o) 

iiear Hc 2 are w^in the observable range. We 
expect to see the enhanced behavior towards Hc 2 : a 
similar behavior already observed in CeCoIns [56] which 
is known to be a typical superconductor with strong PPE. 

(B) To determine the detailed gap structure on the a, 

(3 and 7 bands the field angle resolved specific heat ex¬ 
periment is decisive. The existing data m at low fields 
are only down to 100 mK which was concluded to have 
a dxy like nodal structure, judging from the existing four 
fold oscillation pattern. We expect that the sign of the 
four fold oscillation in C/T might change because the 7 
band responsible for this oscillation at low fields and low 
temperatures should have a like nodal structure. 

Such a sign change of the oscillation has been observed 
in CeCoIns m- This dx 2 -y 2 like nodal structure on the 
7 band is fully consistent with other experiments. In 
particular, the square vortex lattice pattern oriented to 
(110) direction is observed for B || c in the SANS exper¬ 
iments (58). Note in passing that the gap structures of (3 
and a bands are relatively isotropic [65]. This might be 
consistent with the c-axis tunneling data m that probes 
selectively the (3 band with the least Fermi velocity along 
this direction and shows a full gap. 

(C) Finally according to the present theory with PPE, 
Sr 2 Ru 04 is quite likely to exhibit the Eulde-Eerrell- 
Larkin-Ovchinnikov (EELO) phase in low T and high H 
regions just below Hc 2 - One of the best way to detect 
it is to measure by NMR, which is enhanced in the 
EELO phase. This method is successfully applied to the 
organic superconductor /i:-(BEDT-TTE) 2 Cu(NCS )2 [7^ . 

In summary, we investigate the mixed state properties 
of a uniaxially anisotropic superconductor with the Pauli 
paramagnetic effect on the basis of microscopic Eilen- 
berger theory in the clean limit, assuming a single band 
model. By these studies, we discussed the field orienta¬ 
tion dependence near B || a 6 , and tried to explain cu¬ 
rious behaviors in Sr 2 Ru 04 , focusing on contributions 
by the Pauli-paramagnetic effect. In the study of the 
longitudinal and transverse components of the flux line 
lattice form factors, effects of Hc 2 ,ab suppression due to 
the paramagnetic pair breaking is important to quanti¬ 
tatively explain the intensity of the spin-flip SANS ex¬ 
periment observed in Sr 2 Ru 04 [4]. In the magnetiza¬ 
tion curve and field-dependent Sommerfeld coefficient, 
the jumps at the first order Hc2,ab transition mainly come 
from the contribution of the paramagnetic susceptibility. 
Erom the study on the field orientation dependence of 
torque curves and Hc 2 { 0 ), the intrinsic anisotropy ra¬ 
tio between c and ab directions is F ^60, suggesting 
the main contribution of the /3-band in the superconduc¬ 
tivity of Sr 2 Ru 04 . These consistent behaviors between 
experimental observation and the theoretical calculation 
indicate that the Pauli paramagnetic effect plays impor¬ 
tant roles to understand curious behaviors at high fields 
for B II ab in Sr 2 Ru 04 . This suggests that the pairing 
symmetry is either a spin-singlet or a spin-triplet pair¬ 
ing with the d-vector components in the ab plane. We 


expect further experimental and theoretical future stud¬ 
ies to confirm the mechanism of high field behaviors for 
B II a 6 in Sr 2 Ru 04 . 
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Appendix A 

According to Ref. |66l an analytic expression for the in¬ 
dependence of the upper critical field Hc 2 (z^) was derived 
by the fitting to the numerical solutions of Eilenberger 
equation under PPE at low temperatures. The expres¬ 
sion for Hc 2 { 0 ,in) is given by 

Ho2{e,fi) _ 1 

yi + 2.4/i(0)2 

for each field orientation d in a uniaxial superconduc¬ 
tor. There id° 2 ^(< 9 ) = idc 2 (< 9 ,/i = 0) = is the 

orbital limited upper critical field without PPE. The ef¬ 
fective paramagnetic parameter fl{0) at each 0 depends 
on anisotropy T{0) in Eq. (§ as 

m = (A2) 


since /i(d) oc /Hp. The Pauli paramagnetic criti¬ 

cal field Hp is proportional to the pair p oten tial A at a 
zero field. By combining Eqs. (lAll) and (A 2 ), we obtain 


ide2(d,/i) 


1 


\/< 


cos^ 0 -h F”^ sin^ 0 


2. V 


(A3) 


This gives the d-dependent Hc2{0) in Figs. mib) 
and [TT]^c). When we eva luate the T-dependence of 
Hc2{0^ 11 )/in Eq. (A3), = ii{T = 0) is modified 

to the T-dependent paramagnetic parameter /i(T) given 
by 


/i(T) = II 


A(Q) 

A(T) (0) • 


(A4) 


The T-dependent orbital limited given by 

the Werthamer-Helfand-Hohenberg (WHH) formula [75] 
or the solution of the Eilenberger equation in the clean 
limit. The T-dependent order parameter A(T) follows 
the BCS form. The T-dependence of Hc 2 anisotropy 
Hc 2 ah{T)/H c 2 ,c{T) when 0 = 90° is obtained by Eqs. 
(A3) and ( |A4| ), which is displayed in Eig. [TT]^d). 
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